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The extrapolation of self— consistent nuclear mean-field models to the region of 
superheavy elements is discussed with emphasis on the extrapolating power of the 
models. The predictions of modern mean-field models are confronted with recent 
experimental data. It is shown that a final conclusion about the location of the 
expected island of spherical doubly-magic superheavy nuclei cannot be drawn on 
the basis of the available data. 



1 Introduction 

The last decade brought with the synthesis of superheavy nuclei with Z = 110-112 
at GSI (Darmstadt) and JINR (Dubna) a renewal of the interest in the properties 
of superheavy nucleiEm. These are by definition those nuclei with Z > 100 which 
have a negligible liquid-drop fission barrier and are stabilized by quantal shell ef- 
fects only. The ultimate goal is to reach an expected "island of stability" located 
around the nexLsphcrical doubly-magic nucleus which was predicted to be f||ll4 
thirty years agcOaa. Recent experiments performed at JINR DubnaQ give evidence 
for the synthesis of the neutron-rich nuclides 283 112, 287 - 289 H4 ; 292 H6, while at 
Berkeley three a-decay chains attributed to the even heavier 293 118 were observed! 
The measured a-decay. ciiajns of these new nuclides turn out to be consistent with 
theoretical predictionsQ'Bll3. While earlier superheavy nuclei could be unambigu- 
ously identified by their a-decay chains leading to already known nuclei, the decay 
chains of the new-found superheavy nuclei cannot be linked to any known nuclides. 
Tji^ new discoveries still have to be viewed carefully, see e.g. the critical discussion 

The recent experimental developments are accompanied by a significant progress 
in the theoretical modeling of superheavy nuclei by means of nuclear mean-field 
(MF) models. MF models can be divided into two major classes, (i) self-consistent 
(SC) ones where the single-particle wave functions are calculated from an aver- 
age nuclear potential which in turn depends on these wave functions and (ii) 
macroscopic-microscopic (MM) models which are composed by a generalized liquid- 
drop model that governs the bulk properties and a single-particle potential from 
which the shell correction is derived. 

Recent theoretical work reveals systematic differences among the predictions 
of the models. While modern refined MM models confirm the older prediction of 
f |f 1 14 for the next spherical doubly-magic nucleus, nearly all SC models shift that 
property to higher charge numbers, depending on the class of SC models to either 
2 y2l20 or 1I4I26. The reasons for these conflicting predictions and their implications 
are the topic of this contribution. 
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2 Models 



Most SC mean-field models can be viewed as energy density theories in the spirit of 
the Hohenbcrg Kohn Sham (HKS) approachEBE^ originally introduced for many- 
electron systems which is nowadays a standard tool successfully applied in atomk^ 
molecular, cluster, and solid-state physics. Starting point is the existence theoremEJ 
for a unique energy functional £ depending on all local densities and currents that 
can be constructed from the general single-particle density matrix 

p = p{r, a, t- r', o\ t') =^ v l ^( r '> M*, <r, *) (1) 

k 

which gives the exact ground-state energy of a system of identical Fermions when 
£ is calculated for the exact ground-state density p. r, a, and t are the spatial, 
spin, and isospin coordinates of the wave functions ipk- The HKS approach maps 
the nuclear many-body problem for the "real" highly correlated many-body wave 
function onto a system of independent particles in effective, so-called Kohn-Sham 
orbitals ip).. The equations of motion of the ipk are derived from a variational 
principle 

S£ = h(r, a,t)ipk(r, a,t) = e k ip k (r, a,t) (2) 

where the single-particle Hamiltonian h is the sum of the kinetic term t and the 
self-consistent potential T that is calculated from the actual density matrix 

h = ^ = i + r[p} . (3) 
dp 

The existence theorem for the energy functional, however, makes no statement 
about the actual structure of the effective interaction. Guided by symmetry prin- 
ciples and phenomcnological knowledge about nuclei, the aim is to find the most 
simple energy functional which incorporates all relevant physics and to adjust its 
parameters uniquely to a selected set of nuclear key data. The two most widely used 
SC models are the (non-relativistic) Skyrme-Hartree-Fock (SHF) model and the 
relativistic mean-field (RMF) model. The SHF energy functional contains all bilin- 
ear combinations of local densities that are invariant under rotational, translational, 
parity, and time-reversal transformations up to second order in the derivatives plus 
a simple density dependence^. The standard RMF energy functional assumes the 
nucleus to be a system of Dirac nucleons interacting via scalar and vector fields usu- 
ally associated with er, u, and p mesons, again plus a simple density dependencolj. 
Pairing is treated in both SHF and RMF using the same (non=relativistic) local 
pairing energy functional corresponding to a delta pairing forcal3. 

MM models can be .motivated as an approximation to SC models by means of 
the Strutinsky theorcmHEj. The binding energy is separated into a large average 
part E depending smoothly on N and Z and a small shell correction -E s heii that 
describes local fluctuations of the binding energy caused by variations of the density 
of single-particle levels around the Fermi surface 

E(Z,N) = E(Z,N)+E shcll (Z,N) (4) 
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For superheavy nuclei E is of the order —2000 MeV while S s hell fluctuates in the 
range —15 MeV < -E s hcii < +15 MeV. In MM models the self-consistent coupling 
of shell structure and bulk properties is replaced by (independent but similar) 
parameterizations of the the density distribution and the single-particle potentials 
with N and Z . Modern MM models combine a finite-range liquid-drop (YPE) 
or droplet (FRDM) model for E and a phenomenological single-particle model 
based on either the Woods-Saxon (WS) or the Folded- Yukawa (FY) potential for 
the calculation of E^heii- The two most widely used models are the FRDM+FYE3 
and the YPE+WS modelsEJ. It is to be noted that loosely-bound systems like 
superheavy nuclei where the Fermi energy is close to the continuum require a more 
careful treatment of unbound states thap_done in the standard approach used in 
large-scale calculations with MM modelsfea. 



3 Nuclear exotica in superheavy nuclei 



In superheavy nuclei the repulsive Coulomb interaction is not counteracted by the 
surface tension which leads to a vanishing liquid-drop fission barrier. Therefore 
the Coulomb field cannot be treated as a small perturbation atop the nuclear mean 
field, it pushes the protons to the nuclear surface which is counteracted by the 
symmetry energy restoring a similar distribution of protons and neutrons and the 
density dependence of the effective interaction restoring the saturation density of 
nuclear matter (for given asymmetry). With that superheavy nuclei probe the 
balance of bulk properties of effective interactions. At the same time the large 
density of single-particle states provides a sensitive probe for even subtle details 
of the single-particle structure. Again the Coulomb potential induces significant 
changes in the proton shell structure, compared to lighter nuclei single-particle 
states with large angular momentum j are lowered compared to small-j statesEj. 

Up to now most of the understanding of phenomena in superheavy nuclei was 
obtained on the basis of MM models. Recent work employing SC models, however, 
predicts some new phenomena which cannot be consistently described by (current! 
MM models. Examples are exotic radial density distributions like "semi-bubbles"cJ 
(which are essential for the appearance-of spherical shell closures at Z = 120 and 
N = 172 predicted by some SC modelst-3) and the variation of the surface diffusc- 
ness in superheavy nuclei. It is well-known that the appearance of a proton shell 
at Z = 126 requires a larger surface diffusencss than assumed in the standard pa- 
rameterizations of the FY and WS single-particle potentialsEalHa. Standard MM 
models use a parameterization of the radial shape of the density distribution and 
single-particle potentials that has no free parameters to minimize the binding en- 
ergy and therefore are not flexible enough to incorporate both of these effects. It 
can be expected that these models loose their validity for nuclei with large charge 
number. SC models have to be preferred when describing these phenomena as they 
make no assumption at all on the profile of the density distributions of protons 
and neutrons. Besides the limitations of the currently used parameterizations the 
deeper reason why these phenomena cannot be easily described by MM models is 
that they are caused by the coupling of single-particle degrees of freedom and bulk 
properties which becomes more pronounced for loosely bound systems. 
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4 Fits and Parameterizations 



Although they might not contain all physics relevant for superheavy nuclei beyond 
the known region current MM models give a better description of binding energies 
than the best SC modelsEj. This is not too surprising as the aims and the fit 
strategies of MM and SC models are very different: MM models are optimized 
for the description of masses by adjusting the parameters of the macroscopic part 
E to all known masses. The density distribution entering the macroscopic part 
of the model is parameterized as a function of N and Z to reproduce the global 
systematics of radii, the microscopic potential is adjusted to reproduce systematics 
of single-particle spectra throughout the chart of nuclei. On the other hand, SC 
models are designed to give a consistent description of nuclear ground states, single- 
particle spectra, collective excitations like giant resonances and rotational bands, 
and large amplitude collective motion within the same model. The parameters of 
the interaction are usually adjusted to masses and data on the charge distribution 
of selected spherical nuclei, and in some cases to selected data on infinite nuclear 
matter (INM) to compensate for the small number of sample points in the fit. As the 
description of single-particle spectra and bulk properties cannot be separated the 
interactions need not to be adjusted to spectral data with the exception of the spin- 
orbit (£ ■ s) interaction. As a purely relativistic effect it is naturally incorporated 
in the RMFO which reproduces data on_£ • s splittings without being adjusted to 
any data on single-particle spectra at alllLZl. In non-relativistic models - either MM 
or SC ones - the I ■ s force has to be put in by hand and needs to be adjusted to 
spectral data. 

A detail to be kept in mind when comparing masses from MM and SC models 
is that in the fit of MM models the absolute error p£ £ is minimized, while for SC 
models usually the relative error of £ is minimizedEH which allows for much larger 
absolute errors in heavy nuclei, up to 5.5 MeV in 208 Pb corresponding to 0.35%. 

Besides the technical reason that everything else than a fit to a small sample 
of spherical nuclei is too time-consuming there are also physics reasons for the 
usual small sample of fit nuclei. The many-body wave function $ of mean-field 
models breaks symmetries which are obeyed by the effective energy functional and 
the "exact" wave function. An instructive example for this "symmetry dilemma" 
is violation of translational symmetry which is unavoidable as the center-of-mass 
(cm.) of the nucleus is localized by the mean-field potential. Although $ has 
vanishing total momentum (P) = it is not an eigenstate of the momentum oper- 
ator, (P 2 ) ^ 0. This means that $ is not the pure ground state but contains an 
admixture of excited states with finite momentum which decreases the calculated 
mass. Rigorous restoration of the broken symmetry by means of projection is too 
costly to be used in large-scale calculations. The method of choice is to estimate 
the contribution from excited states to the calculated energy and subtract this cm. 
correction i? c . m . to obtain the binding energy. An unexpected side-effect is that 
the nuclear matter properties of effective interactions depend on technical details 
of the actual cm. correction performed during the fit as the effective interaction 
has to compensate for the difference between approximations and the exact value of 
E c . m .. Often-used simple schemes for cm. correction lead accidentally to too large 
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Force 
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&surf 




[MeV] 


[MeV] 


[MeV] 


SkM* 


-15.8 


30.0 


17.6 


SkP 


-15.9 


30.0 


18.0 


SkI3 


-16.0 


34.8 


17.5 


SkM 


-15.9 


29.5 


17.3 


SLy6 


-15.9 


32.0 


17.4 



Force 


Ovol 


G-sym 


Q-surf 




[MeV] 


[MeV] 


[MeV] 


NL-Z 


-16.2 


41.7 


17.7 


NL-Z2 


-16.1 


39.0 




NL3 


-16.2 


37.4 


18.5 



Table 1. Compilation of nuclear matter properties for a number of typical parameter sets. 
SkM*-SLy6 are Skyrme forces, and NL-Z, NL-Z2, and NL3 RMF forces. a vo i denotes the volume 
coefficient or energy per nucleon, o sym the (volume) asymmetry coefficient, and a sur f the surface 
coefficient. Empirical values for the volume cpsfficients derived from the liquid-drop model are 
«vol = -16.0 ± 0.2 and flsym — 32.5 i 0.5. Seet_j for more details. 



a surface coefficientES which has visible impact on the extrapolation of the models 
to heavy systems and large deformation. A similar influence on the results of a fit 
can be expected for other corrections for spurious motions, namely the rotational 
and vibrational correction. Luckily their contribution to the binding energi^can be 
suppressed by choosing spherical nuclei with stiff potential energy surfacesEH which 
explains the usual small sample of fit nuclei. 

There are numerous parameterizations of the SHF and RMF models to be found 
in the_Jiterature. Resplts discussed here_are obtained with the SHF interactions 
SkM*H, Skip, SLy6B, SkI3 and SkI4B, and the RMF forces NL3B, NL-Zo, 
and NL-Z2tLZl. SkP uses effective mass m* /m = 1 and is designed to describe both 
mean-field and pairing interaction. The other Skyrme forces all have smaller effec- 
tive masses around m*/m ~ 0.75. SkM* was first to deliver acceptable incompress- 
ibility and fission properties. SLy6 stems from an attempt to cover properties of 
pure neutron matter together with normal nuclear ground-state properties, while 
SkI3 and SkM stem from a recent fit including data from exotic nuclei and use a 
variant of the Skyrme parameterization where the I ■ s force is complemented by an 
explicit isovector degree-of-freedomo. They are designed to overcome the different 
isovector trends oil - s coupling between conventional Skyrme forces and the RMF. 
The RMF parameterizations NL-Z, NL-Z2, and NL3 use the standard non-linear 
ansatz for the RMF model. NL-Z aims at a best fit to nuclear ground-state prop- 
erties. NL-Z2 matches exactly the same enlarged set of data like the Sklx forces. 
NL3 results from a fit including neutron rms radii and nuclear matter data. 



5 Predictive Power 

Table | summarizes those nuclear matter properties that can be directly linked to 
the leading terms of the liquid-drop model 

-Bldm = a vo i A + a sym I 2 A + a surf A 2/3 + ... (5) 

for some typical SHF and RMF forces. Most non-relativistic (relativistic) interac- 
tions agree among each other in the values for a vo \ and a sym , but relativistic and 
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Figure 1. Relative error SE = (S ca i c — i? C xpt ) /-^oxpt (in %) of the binding energy for the chain 
of Pb isotopes and the chain of N = 126 isotones. Negative values denote under-bound, positive 
values over— bound nuclei. |Tlji | u h orizontal lines at ±0.35% indicate the relative error in binding 
energy allowed for good fitsEjl23. 



non-relativistic interactions differ significantly. While the Skyrme force values are 
close to the empirical ones, the R,MF forces give larger d V oi ^rid tJ S y m . work much 
better in that respect. It is important to note, Not all values for Skyrme interac- 
tions are predictions, for the SLy6, SkP, and SkM* forces nuclear matter properties 
were used as input data during their fit. As already mentioned above the predicted 
value for a sur f is accidentally correlated to the scheme for cm. correction used in 
the fit of the interactionEa, forces with simple cm. correction like SkP and NL3 
have systematically larger a sur f than the others. 

A first impression of the quality of mean-field models for finite nuclei is provided 
by the systematics of the relative error of binding energies, see Figs. [j] and || for 
typical results. For the chain of spherical Pb isotopes the results remain essentially 
within the bounds of 0.35% allowed in the fits, while the spherical TV = 126 isotones 
reveal already some unresolved trends in the RMF interactions. The slope of 5E 
corresponds to an error of the two-nucleon separation energies. Unlike modern 
SHF interactions which have four or even five independent isovector terms, the 
RMF has one isovector term only which seems to be insufficient. The small change 
in the isovector coupling between NL-Z and NL-Z2 improves visibly the trends 
with N but at the same time worsens the trends with Z. But one has to be careful 
as both I and A change along isotopic and isotonic chains which mixes trends in 
the isovector and isoscalar channel. These can be separated plotting SE against 
I and A, see Fig. |2| with data for superheavy nuclei. It can be expected that in 
this extrapolation of the models the SE spread more than in the case of the lighter 
nuclei. It is gratifying, however, to see that most of the interactions stay within 
or at least close to the bounds of 0.35% error. When plotted versus A one obtains 
essentially flat curves for all SHF forces, while there is a small but visible slope 
for the RMF forces that points at an unresolved isoscalar trend. When plotted 
versus / all forces (perhaps with the exception of SkP) show slopes which point 
at unresolved trends in the isovector channel, even for modern forces like NL3 or 
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Figure 2. Relative error the binding energy for the heaviest deformed nuclei where the mass is 
known plotted for constant I = N — Z against A = N + Z (left panel) and vice versa (right panel) 
to separate trends in the isosc aU i j nd isovector channels of the effective interactions. Symbols as 
in Fig. hj. Data are taken fromLJCJ. 



SLy6 fitted with bias on good isovector properties. One has to be careful with the 
interpretation of Fig. ^| since the range of known masses in / is rather small and 
wrong trends in the "macroscopic part" of the SC models might interfere with local 
flaws in the microscopic part. An example for that is the peak in the 8E around 
208 Pb in the left panel of Fig. [I] which hints at an unresolved shell effect. It is to be 
noted that missing shell effects might also correspond to correlations beyond the 
mean-field level, see the example of cm. correlations discussed inn3. 

As there are no experimental data on spectral properties of superheavy nuclei, 
the predictive power of the models for single-particle spectra has to be examined 
looking at lighter nuclei, see Fig. ||. None of the existing parameterizations of 
SC or MM models is able to give a proper description of the level ordering and 
the spin-orbit splitting in heavy nuclei such as 132 Sn or 208 Pb. Levels with large 
angular momentum are usually pushed up too far in the single-particle spectrum 
and all non-relativistic models (either SC and MM ones) show a wrong trend of the 
spin-orbit splitting with A and therefore usually overestimate the proton spin-orbit 
splitting in heavy nuclei. This is devastating for interactions where the spin-orbit 
coupling constant is adjusted to data for 16 0. Including data on heavy nuclei in 
the fit gives better overall agreement but cannot remove the overall wrong trend 
as can be seen from the example of SkP. It is surprising that Skyrme forces with 
an additional degree of freedom in the spin-orbit interaction (SkI3, SkI4) perform 
worse in that respect than standard forces (SkP, SLy6). The overall performance 
of the RMF interactions is much better. 

As most mean-field models predict the same level-ordering in the superheavy 
region, already slight changes in the relative distances among the interactions lead 
to different magic numbers, see Fig. |^ for typical results. Large spin-orbit splitting 
as in case of the FY and SkI4 models favors Z — 114, but as these interactions 
overestimate the spin-orbit splitting of proton states in heavy nuclei this prediction 
is very doubtful. Non-relativistic models in general prefer N — 184 for neutrons, 
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Figure 3. Single— particle spectra of the protons in 208 Pb (left panel) and relative error of the 
spin-orbit splittings Se ls = (ei SiCalc - ei s ,cx P t)Als. cxpt ( m %) °f proton states close to the Fermi 
surface in ls O, 132 Sn and 208 Pb for the interactions as indicated. Negative errors denote calculated 
values which are too small. FY denates the folded- Yukawa single— particle potential widely used 
in MM models. Data are taken fromO. 



while small spin-orbit splitting in connection with "semi-bubble" shapes leads to 
Z = 120 and N = 172 as it happens for SkI3 and the relativistic forces NL3, NL-Z, 
and NL-Z2. SC interactions with large effective mass close to m* jm w 1.0 like SkP 
give an average level density so large that most of the shell effects are washed out 
and the proton shell is shifted to Z = 126. A common feature of spherical shells 
as predicted by SC models is their strong nucleon-number dependencetZfED: While 
the spectra for the non-SC FY model remain basically unchanged, the spectra from 
SC models change dramatically going from Iff 114 to ffflSO, the spherical shells at 




Figure 4. Single— particle spectra of the protons in f§f 114 and 2 7?,120 at spherical shape for the 
interactions as indicated. Note that in spite of the|-dj>»ble shell closure the ground state of 120 
is actually deformed for most Skyrme interactions£3o, while the strongr^V = 184 shell stabilizes 
for most interactions the spherical shape of 2 g|ll4. Data are taken fromtj. 
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Figure 5. Proton shell correction extracted from calculated self— consistent binding energies of 
the N = 172 and N = 184 isotones (left and middle panel) and neutron shell correction for the 
Z = 120 isotopes. The same energy scale is used for all panels. Note that for the majority 
of Skyrme interactions all N = 172 isotones are predicted to be deformed. The inserts show 
the density profile for 292 120 (left panel) and 304 120 (middle panel) calculated with NL-Z2. 
They demonstrate the strong neutron-number dependence of the "semi-bubble" i-shapes which 
are responsible for the magic numbers Z = 120 and N = 172. Data are taken from! 22 !. 



Z = 120, Z = 126 and N — 172 are all restricted to a narrow range of N and Z, a 
magic N — 184 often excludes a shell closure at Z — 120. 

The appearance of gaps in the single-particle spectra alone is not sufficient to 
stabilize a superheavy nucleus. Much more important is the extra binding and its 
strong shape dependence obtained from a smaller than average level density. The 
shell correction can be viewed as the natural measure for this "shell effect" and 
provides a powerful tool to analyze also results obtained in fully self-consistent 
calculations, see Fig. [|. It has to be emphasized that the shell correction is not 
equivalent to the gap in the single-particle spectrum and that a quantitative com- 
parison of the two quantities cannot be made. For N — 172 most of the Skyrme 
forces give a magic Z = 120, while for N = 184 the minimum is shifted to Z = 124- 
126 in all cases. RMF forces give consistently different results, the strongest shell 
effect is at Z = 120 independent on the neutron number. Most interactions pre- 
dict also a broad region of negative neutron shell correction. Again a systematic 
difference between the models: while the RMF points at TV = 172 as the domi- 
nant neutron shell, all SHF forces predict N — 184. Note that shell corrections 
for neutrons are much larger than for protons. This partly explains the finding 
that spherical ground states of superheavy .nuclei are usually correlated with magic 
neutron numbers, not proton numbersE3E3il3. All this suggests that the "island of 
stability" is not so much coupled to particular shell closures but a region of nuclei 
with low level density which is also found in MM modelsEd. Fig. [| demonstrates 
also that it is nearly impossible to discriminate between models giving different 
shell closures by looking at systematics of binding energies or Q a values alone. Al- 
though the proton shell correction might be peaked at different points, its variation 
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Figure 6. Macroscopic energy E extracted from calculated self-consistent binding energies (left 
panel) and macroscopic energy Sldm from the LDM expansion (H) of the binding energy (right 
panel) for the chain of N = 184 isotones. The phenomcnological macroscopic energy from the 
YPE model is given for comparison. Note that the scales of the two panels are very different. To 
illustrate *he Z dependence all energies are normalized to the YPE value at Z = 100. Data taken 
from Ref.B. 



is often too small to be visible within the uncertainty of SC models. 

The shell correction does not only extract the "shell effect" from the self- 
consistent binding energies, at the same time one obtains the macroscopic part 
E of the binding energy. A typical example is given in the left panel of Fig. [| 
Owing to its optimization as a mass formula and its fit which includes also super- 
heavy nuclei the YPE model can be expected to give a very good description of the 
macroscopic energy in this region and therefore can serve as reference. The curves 
for the various interactions show considerable splitting which reflects of course the 
findings from Figs. |l| and 0. This also confirms again the finding discussed above 
that the predictive power of an effective interaction for binding energy systematics 
is fairly independent of its predictive power for shell effects. On one hand SLy4 and 
NL3 give very similar values for E as the YPE model (but remember the small er- 
rors in their trends found in Fig. ^) although they predict different magic numbers, 
on the other hand those forces that give the best description of binding energies 
independently on their disagreement on the next shell closure (SkI3, SkI4, SkP, 
NL-Z2) give all very similar E but differ now significantly from the YPE model. 

The right panel of Fig. ^ shows the macroscopic energy -Eldm from the LDM 
expansion of the binding energy (|5|) using the values given in Table |l|. Surprisingly 
E and -Eldm differ on the order of 100 MeV. Only for -Eldm the ordering of the 
interactions is according to their value for a sym as one naively expects. This demon- 
strates that the leptodermous expansion with nuclear matter parameters does not 
work even for superheavy nuclei, finite-size effects are still important which also 
means that the use of nuclear matter parameters as pseudo-observables in the fit 
of SC models might be dangerous. 

Most of the recent new data on superheavy nuclei are Q a values. Their sys- 
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Figure 7. Q a values of even-even nuclei calculated with NL-Z2 (open circles) compared with 
experimental values for even-even nuclei (filled diamonds) and odd— N isotopes (open diamonds) 
of the heaviest even Z elements. The data for odd-TV nuclei have to be handled carefully, some of 
these might correspond to transitions involving excited states, and due to blocking effects of the 
ground— state-to— ground— state values might differ on the order of 500 keV from the systematics 
of Qa for even-even nuclei. Including data for nuclei wift Z > 116 leads to overlapping curves, 
therefore those are omitted in the plot. Data taken fromPT 



tematics reflect all properties discussed above. As a-decay chains have / =const. 
the isoscalar channel of E mainly determines the overall slope of the Q a , while the 
isovector channel of E mainly shifts the whole curves around. Shell effects bend the 
curves locally, leading to kinks and peaks. A model has to give a perfect description 
of all these properties to reproduce experimental data throughout the superheavy 
region. The quality of NL-Z2 for the Q a is shown in Fig. [?]. The overall description 
of the data is very good with the exception of some nuclei around Z = 104 where it 
overestimates a deformed shell closure while the deformed N = 152 shell is shifted 
to N = 150. The latter is a problem from which virtually all SC models sufferBEl. 

Fig. H compares some very recent data with calculations. In view of the un- 
certainties, the SC SLy4 and NL-Z2 give a very good description of the data for 
the decay chain of iHll2 (NL-Z2 agrees for Z > 104 only as can be expected 
from Fig. [?]) and reproduces the N = 162 shell effect which cannot be seen in the 
FRDM+FY predictions. While all models give similar predictions for this well- 
established chain, the spread among the models is much larger for the new chains. 
All models with the exception of YPE+WS show spherical or deformed shells which 
cannot be seen in the data. 

Comparing predictions with the recent data for the even-even ^§116 decay 
chain (which still have to be viewed as preliminary), see the right panel of Fig. ||, 
it is most interesting that the data agree with calculated values from interactions 
which give different predictions for the spherical magic numbers, i.e. SkI4 (Z = 114, 
N = 184), SLy6 (Z = 126, N = 184), and NL3 (Z = 120, N = 172). All other 
interactions show wrong overall trends of the Q a or pronounced deformed shells 
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Figure 8. Left panel: Comparison of experimental and calculated Q a values for the decay chains 
of ^^112, ^y|ll4, and ^118, in the latter two cases following the mass and charge assignment 
of the experimental groups. In the 1I5II2 chain two distinct branches leading through different 
states of the intermediate nuclei are known. The calculated values from NL-Z2 and SLy4 connect 
the lowest states with positive parity in all cases (in the new chains only ^|114 and ig^HO are 
predicted by SLy4 and NL— Z2 to have ground states with negative parity), while the FRDM+FY 
and YPE+WS data are ground state to ground state values. Right panel: Q a values of nuclei in 
the decay— chain jggl20 — » 1 IS — * . . . jggHs as predicted by the mean— field models as indicated 
compared with preliminary experimental data. 



in disagreement with the data or even both. The large difference between NL-Z 
and NL-Z2 is caused by their difference in symmetry energy although both forces 
predict the same shell structure. Once again all this demonstrates that predictions 
for spherical shell closures and binding energy systematics are fairly independent. 

One of the reasons for the poor description of the new data is that - independent 
on the actual location of the shell closures - these nuclei are located in a region 
of transitional nuclei with very soft potential energy surfaces which amplify small 
differences in the shell structure. At the same time this adds a large uncertainty 
to the predictions as correlations beyond the mean-field level give a non-negligible 
contribution to binding energy differences which washes out shell effects visible in 
the results from mean-field calculations^!!. 



6 Conclusions 



Combining the findings for average trends of binding energies and single-particle 
spectra it has to be said that none of the current models is able to describe all 
available data. Errors in the macroscopic part of the models have to be distin- 
guished from errors in the shell structure. A good description of binding energy 
trends is not neccesarily an indicator for the predictive power concerning shell ef- 
fects. Looking at "macroscopic" observables, a single parameterization of a model 
is not representative for the model while observables sensitive to shell effects re- 
veal systematic differences between SHF and RMF. The good overall description 
of spin-orbit splittings by the RMF gives some preferrence for its predictions for 
magic numbers, but the undoubtedly missing isovector degrees of freedom in the 
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RMF once included might also feed back to its predictions for superheavy shell 
closures. More research in that direction has to be done. 
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